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ABSTRACT 

The core of the nearby galaxy NGC 660 has recently undergone a spectacular radio 
outburst; using a combination of archival radio and Chandra X-ray data, together with 
new observations, the nature of this event is investigated. Radio observations made 
using e-MERLIN in mid-2013 show a new compact and extremely bright continuum 
source at the centre of the galaxy. High angular resolution observations carried out 
with the European VLBI Network show an obvious jet-like feature to the north east 
and evidence of a weak extension to the west, possibly a counter-jet. We also examine 
high angular resolution Hi spectra of these new sources, and the radio spectral energy 
distribution using the new wide-band capabilities of e-MERLIN. We compare the 
properties of the new object with possible explanations, concluding that we are seeing 
a period of new AGN activity in the core of this polar ring galaxy. 

Key words: galaxies: individual: NGC 660 - techniques: interferometry - techniques: 
high angular resolution - radio continuum: galaxies 


1 A POLAR RING GALAXY 

The nearby galaxy NGC 660 (~13Mpc; ISpringob et al.l 
I2OO9II exhibits a polar ring morphology and a LINER- 
type nucleus (low - ionisation nuc l ear em ission-line region; 
iNagar et al.lI2OO2I I . ICombes et al.l lll992tl used observations 
of CO emission to infer a total H2 mass in the main disk 
(radius 4kpc) of 3.7 xlO® M©, and in the polar ring (radius 
12kpc) of 10® Mq. They construct a model to account for 
both the CO and Hi properties of the system, concluding 
that the ring is highly inclined and highly warped. The stel¬ 
lar popul ation in the ring is es timated to be a few billion 
years old (Ivan Driel et al.l 19951 ) . 

The large-sca le Hi and OH properties were probed by 
iBaan et al.l lll992l l who saw a large velocity gradient from 
the disk in absorption, as well as absorption from the warped 
outer disk at the systemic velocity. The same authors also 
detected a molecular outflow in the core, thoug ht to be due 
to a p erturbed spiral structure in the inner disk. iBaan et al.l 
lll986ri also searched several galaxies for an H 2CO signal, dis¬ 
coveri ng weak absorption in NGC 660, while IPhookun et al.l 
lll998ll looked at the H92a radio recombination line and 
found the dynamical centre is offset from the continuum 
peak. More recently, H2CO and NH3 absorption have both 
been detected against the central region of the galaxy but 
with very different widths, suggesting that the two lines 


are pr obing different parts of the galaxy dMangum et al.l 
l2013all^ . 

The Hi and CO distributions have been examined at 
low spatial re solution, as wel l as th e optical and infrared 
properties, by Ivan Driel et ^ (Il995l l who found a flat disk 
rotation curve, a high average dust temperature similar to 
that of M82, but a low star-formation efficiency more similar 
to that found for non-starburst spirals. The main gas-rich 
disk cont ains compact (on scales of 0'.'2) radio continuum 
emission llCarral et al.|[l9^ 1. There is weak extended radio 
emission from the centre of the galaxy associated with the 
main disk, and observations with the VLA in C-array at 
8.4 GHz show many co mpact objects within the disk of NGC 
660 dFilho et al.ll2002lb 

Low angular resolution archive radio observations show 
NGC 660 to be a fairly unremarkable starburst galaxy with 
no obvious active nucleus (AGN), and any radio emission 
from pre-existing compact so u rces appears to be purely 
starburst-related. iDudik et al.l d2005l l investigate a sample 
of IR-bright LINERs, including NGC 660. From the stellar 
velocity dispersion, they estimate the black hole mass in 
NGC660 to be logMen = 7.35 Mq, and calculate the Ed¬ 
dington ratio to be Lboi/LBdd = 4.0 x 10“®, typical for LIN¬ 
ERs and indicating a radiatively inefficient accretion flow. 

W hilst several radio s ources were seen in the 15-GHz 
map of ICarral et al.l d 19901 ) (see also Figure [T]) no compact 
emission was seen at the dynamical centre of the galaxy. At 
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hig her resolution, no e mission from the galaxy was detected 
by Ijones et alj (Il98lll in a three-station VLBI experiment 
at 18 cm in 1978. Archival MERLIN data from 1998 (Fig.[T|) 
show no evidence for a compact core at 5 GHz, just a single 
compact (-^SOmas) source with an integrated flux density 
of 2.6 mjy, located 2.5 arcseconds away in the north-eastern 
end of the ring structure. 

At radio wavelengths, the galaxy flared spectacularly 
sometime between 2008 and 2012 and has been reported 
to have developed a new continuum source with a GHz- 

! ieak spectrum and a peak flux density of ~0.5 Jy at 5 GHz 
Minchin et alJl201,'j f. High-resolution imaging has report¬ 
edly shown this increase to be due to a new source in 
the core, the nature of which is currently unclear. Simul¬ 
taneously, OH emission and absorption in both the ground 
state (18-cm) and excited (6-cm) lines has developed, and 
H2CO absorption is seen against the new continuum source 
llMinchin et al.l[2013h . 

Whatever the explanation for the outburst in this 
galaxy, the event is clearly unique. In order to investigate 
its nature we examined archival X-ray and radio data, ob¬ 
served the galaxy with e-MERLIN at both 1.5 and 5 GHz, 
and at 1.4 GHz with the EVN, obtaining both a high an¬ 
gular resolution continuum image to look for evidence of 
small-scale structures, and spectral line data to probe the Hi 
and look for evidence of outflowing gas; this paper describes 
these observations. Section [2] outlines the observational se¬ 
tups used to obtain new data and the archival data which 
has been examined, sections |3] and [4] describe the continuum 
and spectral line results, we explore possible scenarios in[5l 
and summarise our findings in section [S] 


2 OBSERVATIONS 

NGG660 was observed at 1.4 GHz with 11 stations of the 
EVN (Effelsberg, Westerbork, Lovell, Onsala, Medicina, 
Noto, Torun, Svetloe, Zelenchuk, Badary, Urumqi) on 30 
October 2013 (project code EA054). The observation used 
3C84 and 3C120 as bandpass calibrators, and the nearby 
source J0143+1215 as a phase calibrator with a switching 
cycle of seven minutes on NGC660 to three minutes on 
the phase calibrator. The phase calibrator has a position of 
RA= 0l‘'43m3E092221 Dec= -hl2°15'42'.'93343 with po¬ 
sitional uncertainties of 1.02 mas in RA and 0.78 mas in 
Dec. The data were correlated in two passes, one pass of 8 
sub-bands (each with 32 channels, dual polarisation and 2- 
second integrations) for continuum, and one pass with one 
high-frequency resolution sub-band (1024 channels, dual po¬ 
larisation, 2-second integrations) to provide spectral line in¬ 
formation with a velocity resolution of ~3.3 km/s. No serious 
problems were reported during observations, however Noto 
recorded no fringes in the line pass, and several antennas suf¬ 
fered from significant interference. Data reduction proceeded 
according to the EVN User Guide, including: interference re¬ 
jection, fringe fitting, ionospheric corrections, phase calibra¬ 
tion, amplitude calibration and bandpass calibration. Both 
the line and continuum datasets were processed in the same 
way, up to the bandpass stage. Once calibrated, the contin¬ 
uum data were averaged and imaged, whilst each channel 
of the line data was imaged individually to form a cube for 
spectral analysis. Self-calibration was carried out on both 


the phase calibrator and target, initially in phase-only mode, 
then using both amplitude and phase. Maps were made dur¬ 
ing each iteration of self-calibration, and the solutions were 
carefully inspected at each step. 

Lower angular resolution e-MERLIN observations were 
also made in May/June 2013 at both 1.4 and 5 GHz. These 
observations used all six telescopes of the standard array 
(Mk2, Pickmere, Darnhall, Knockin, Defford, Cambridge) 
and used the same phase calibrator as the EVN observa¬ 
tions. For flux and bandpass calibration, 3C286 and OQ208 
were observed using the same configuration. T he data wer e 
processed using vO.7 of the e-MERLIN pipeline llArgdl2015ll . 
with further editing and self-calibration cycles carried out by 
hand. 


2.1 Archival observations 

MERLIN 1.4- and 5-GHz observations from 1998 have been 
extracted from the archive and re-processed, and two Very 
Long Baseline Array (VLBA) observations taken in 2001 
(BF068; 5.0 GHz) and 2010 (BC191; 8.4 GHz) have also been 
examined. The 5-GHz 2001 VLBA data cover a total of 
1^23™ split into two observations, separated by four hours, 
using four 8-MHz sub-bands each split into 16 channels. 
The 8.4-GHz 2010 data consist of two observations of eight 
minutes each, separated by 96 days, each observation was 
made using 8 sub-bands of 16 MHz split into 128 channels 
each. In add i tion, t he VLA A-array 1986 15 GHz data from 
ICarral et al.l (Il990l l has also been extracted and analysed. 

Four separate X-ray observations from the Chandra 
archive were used, taken in January 2001, February 2003 
and, December and November 2012 (ObsIDs 1633, 4010, 
15333 and 15587). The data consist of a total of ~52 ks of ex¬ 
posure time, ~7 ks before the radio outburst and ~45 ks af¬ 
ter. All of these data were obtained using the Advanced CCD 
Imaging Spectrometer (ACIS-S). The data were processed 
using CIAO v4.5, w ith the most recent calibration files 
llFruscione et al.ll2006l ). Spectra covering 0.1-10 keV were ex¬ 
tracted using a circle of radius 6 arcsec centred on the lo¬ 
cation of the radio nucleus, but avoiding the point source 
to the west of the nucleus (see Section Background 

spectra were extracted using an annulus surrounding the 
source extraction region, with an outer radius of 25 arcsec 
and an inner radius of 13 arcsec, avoiding nearby X-ray point 
sources. 

3 CONTINUUM COMPONENTS 

3.1 Radio results: a time-line 

No central compact radio source was detected within 
NGC660 in either archival MERLIN observations from 
1998 (to a 3-a limit of 0.3mJy/bm, Figure [T] left), or 
VLBA observations taken in 2001. The archival VLBA ob¬ 
servations in September 2010 show that a single compact 
source (beam size 2x1 mas) was present at the centre of 
the galaxy at 8.4 GHz by late 2010. The peak flux den¬ 
sity rises from 18.3±0.9 mjy/beam on September 16th to 
36.0±1.8 mjy/beam on December 21st 2010. 

The two e-MERLIN continuum datasets from mid- 
2013 both show a strong compact source at the same po- 
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Figure 1. MERLIN archive data at 4.8GHz from an observation on Jan 9*^ 1998 (left), and VLA A-array 15-GHz data from 1986 
(right; ICarral et al.llT993i . The only radio source visible i n the MERLIN data is a supernova remnant, located towards the north-east 
end of the main disk of the galrtxy. In the VLA data from ICarral et ^ ill99011 the supernova remnant is also present, along with some 
emission from the central star-forming ring. In both images the cross marks the position of the new continuum source reported here. In 
both maps the peak flux is 1.7 mjy/beam; the contours are set at (-1, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10) X 345 ^Jy/beam (30-) in the MERLIN 
map, and (-2, 2, 3, 4, 5, 6, 7, 8, 9, 10) X 180/rJy/beam in the VLA map. 


sition as the 2010 VLBA detection, unresolved in both fre¬ 
quency bands (150mas at 1.4GHz, 50mas at 5GHz). In 
the 2013 e-MERLIN data, the integrated flux density in the 
centre of the 5-GHz band is 364±18mJy. Gomparing this 
with the MERLIN pre-outburst limit results in a factor of 
>1200 change in brightness. At the centre of the e-MERLIN 
band observed at 1.4-GHz, the integrated flux density is 
198±9.9 mjy, giving a fairly flat overall spectral Index of 
a = 0.5, where S oc v°‘. 

The most recent observation, the October 2013 EVN 
data (with an angular resolution of 26x20mas), shows an 
obvious jet-like feature to the north-east of the central 
brightest source (Figure [2l greyscale). When the EVN data 
are imaged with uniform weighting, giving greater resolu¬ 
tion (19x6 mas) at the expense of increased noise, there 
also appears to be a weaker feature to the west (Figure [2l 
contours). The high resolution image shows a compact cen¬ 
tral object with an integrated flux density of 69.1±3.5mJy. 
Despite having greater angular resolution, the 2010 VLBA 
data shows no sign of the structures seen in the EVN data 
taken three years later, only the central brightest source. 
This could be because the object is evolving spatially in 
time, or simply due to the lack of sensitivity and poor u — v 
coverage of the earlier VLBA snapshot observations. 

Between e-MERLIN and EVN scales there is a sig¬ 
nificant decrease in the flux detected. Summing the fitted 
components in the EVN model gives a total flux density 
of ~105±5 mJy, while the e-MERLIN flux density mea¬ 
sured in the sub-band closest to the EVN measurement is 
196±10mJy, almost twice that detected by the EVN. This 
large difference in flux is likely due to a combination of the 
two arrays not being sensitive to the same spatial scales, and 
the source possibly fading as it evolves in the five months 
between the observations. 


Table [T] shows the results of two-dimensional Gaussian 
fits to the EVN continuum images. The two-component fit 
was carried out on the naturally-weighted image, while the 
three-component fit is for the higher resolution uniform- 
weighted image. If component 2 is, as it appears to be, 
moving, then we can estimate limits on its v elocity. Assum¬ 
ing t hat the event occurred between 2008.0 llMinchin et al.l 
l2013ll and 2010.7 (epoch of the first of the VLBA detections), 
then the apparent speed of component 2 is 0 . 6 c ^ i> ^ 1 . 2 c, 
relative to component 1 . 

3.2 Archival X-ray results 

The X-ray image of the nucleus after the radio outburst (0.1- 
lOkeV, Fig.|3)l shows a region of diffuse emission around the 
position of the compact radio core. The emission is extended 
-^8 arcsec to the north-east of the core and arcsec to the 
south-west, coincident with both the orientation of the main 
disk of the galaxy, and the jet-like components seen in the 
high-resolution radio data. The extension seen perpendicular 
to the main disk is not aligned with the axis of the star- 
formation ring and m ay be associated with a starburst wind 
like that seen in M82 (IStrickland et al]| 19971) . There is also a 
peak in the emission at the position of the radio core, which 
is 2-3 times brighter than the surrounding emission, but it is 
unclear whether this is a brighter patch of diffuse emission 
or is actually a peak of emission from the nucleus. 

Unfortunately, the data taken before the outburst are 
not sufficient to make a meaningful image, as the very low 
count rate (1.30 ± 0.10 x 10~^s“^) meant that the exposure 
time was insufficient. The 0.5-10 keV X-ray fluxes (depen¬ 
dent on the model used) were 1.241 q 54 x 10 “^® ergs cm“^ s“^ 
and 1.85lQ jg X 10 “^® ergscm“^s“^ before and after the 
burst, respectively. Assuming a redshift of 0.002842, this 
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Figure 2. Continuum image of the new source in NGC 660 show¬ 
ing the jet-like components. The “low" resolution (natural weight¬ 
ing; 26x20 mas) image is shown in grey scale, with contours from 
the “high" resolution (uniform weighting; 19x6 mas) image over¬ 
laid. The numerals correspond to those of the spatial fit param¬ 
eters listed in Table ^ The data are contoured at (-1, 1, 2, 4, 
8, 16, 32, 64) X 1.736 mjy/beam, the greyscale range is -1.26 to 
66.49 mJy/beam. 



RA 

Dec 

Peak 

Int. 

Tb 


(J2000) 

(J2000) 

(mjy/bm) 

(mjy) 

(xlO'^K) 

1 

02.319 

44.9033 

75.4 ±3.8 

80.0 + 4.0 

>7.2 

2 

02.320 

44.9130 

26.3 + 1.3 

25.6 + 1.3 

>2.3 

1 

02.319 

44.9028 

62.8 + 3.1 

69.1 + 3.5 

>11.6 

2 

02.320 

44.9127 

24.9 + 1.2 

27.7+1.4 

>4.6 

3 

02.319 

44.9068 

10.1 + 0.5 

10.6 + 0.5 

>1.8 


Table 1. Results of 2-component (top) and 3-component (bot¬ 
tom) fits to the naturally-weighted and uniform continuum images 
respectively. Labels correspond to those in Figure [2] Positions are 
relative to 01^43™ 13°38^, and phase referenced to J0143+1215 
(RA = 01^43"^3R092221 Dec = 12°15'42793343; J2000) which 
has a catalogued positional uncertainty of ±1.02 mas in RA and 
0.78 mas in Dec. Brightness temperatures are calculated using 
the beam major axis and are lower limits only. 


corresponds to luminosities of 2 . 20 ^o 50 ^ ergss“^ and 

3.28lo'28 ^ ergss“^ respectively, so potentially brighter 
afterwards. 


4 SPECTRAL FEATURES 
4.1 Radio Hi spectra 

Figure [ 4 ] shows two radio spectra taken from the EVN line 
dataset. The line pass has a higher noise level than that 


Q 



Figure 3. Combined radio (archive 1998 MERLIN 1.4-GHz; 
greyscale) and X-ray (Chandra 0.1-10 keV; contours) image of the 
diffuse emission around the nucleus of NGC 660. The X-ray con¬ 
tours use the combined data from the two observations from 2012, 
after the radio outburst. The contour levels are at (4.5, 12.5, 22.5, 
30.0, 37.5) xl0“® counts s"*^. The greyscale range is from -1.0 
to +4.6 mjy/beam. The X-ray emission is well aligned with the 
radio emission on this scale, which is extended in the direction of 
the main disk of the galaxy. 
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of the continuum dataset, with the result that it was not 
possible to make a reliable spectrum of source 3, the weakest 
of the three sources shown in Figure [2l In Figure IH the solid 
black line is from source 1, while the dashed red line is from 
source 2 (see Figure [2]). It is clear that both radio spectra 
are composed of several narrow spectral components, with 
significant differences between the two. 

The strongest component in the spectrum of t he core is 
very close to the systemic velocity (845±1 kms~^: ILu et ahl 
I 1993 I) . The spectrum from source 2 is also composed of sev¬ 
eral narrow velocity components, but with the features on 
the blue side of the systemic velocity significantly shifted to 
lower velocity, the same sharp peak near 900kms“^, and no 
broad feature near 960 km s”*^. Attempts were made to fit 
the complex line structure with a range of multi-Gaussian 
component models. The best constrained of these models 
comprised of six Gaussian components, the results of which 
are presented in Table[2]and shown in Figure[Sl The reduced 
for this fit is poor (9x10®), but it should be noted that 
excluding any of these components (or adding additional 
Gaussians) made the quality of the fit significantly worse. If 
the new compact radio sources are, as they appear to be, at 
the centre of NGG 660, then each absorption profile is prob¬ 
ing a very narrow line of sight through the near-side of the 
galactic disk, with each of the narrow spectral components 
corresponding to a different physical feature along the line 
of sight. 

The EVN spectrum taken against the brightest source 
(component 1) matches very well with that seen in recent 
Westerbork observations (van Bemmel et al in prep). How¬ 
ever there is less absorption in the red-shifted wing centred 
around 960 km s”*^ than is seen in the spatially-unresolved 
Westerbork spectrum, and no evidence of the blue-shifted 
wing below 850kms“^ seen in lower spatial resolution data. 
This is not surprising, since the EVN spectrum is probing 
a much narrower line of sight through the galaxy than that 
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HI absorption against core and jet in NGC660 



velocity / km/s 



Figure 5. Spectrum against component 1 showing the six Gaus¬ 
sian components from the fit in Table [2] (red dotted lines). The 
composite fitted spectrum is also overlaid (pale dashed line). 


Figure 4. Hi spectrum against components 1 (black, solid) and 
2 (red, dashed) from the October 2013 EVN observation. The 
vertical line marks 845 km s“^ the systemic velocity of the galaxy. 


of the much larger Westerbork beam. Significantly, the deep 
absorption seen on the blue side of the systemic velocity 
in the spectrum taken against source 2 is systematically 
> 5kms“^ (> 4 . 80 ") bluer than that taken against the cen¬ 
tral radio source. It is likely that this is not pur ely due to the 
large-scale disk rotation, as iBaan et al.l ll 19921 ) found a gra- 
dient of 1.58kms“^ pc“^ across the disk in Hi absorption; 
this would result in a difference of just 1.9kms“^ between 
the positions of sources 1 and 2 . 

The EVN Hi spectrum shows that the peak 21-cm opac¬ 
ity, r, is 1.58, f rdv = 66 kms“^ and hence Nhi = 0.182 x 
10^^ (Ts/IOOK) f rdv = 12.1 x 10^^ cm“^ where we assume 
Tg = lOOK. Compared to the results of iGallimore et al.l 
lll999l) on S eyfert and starburs t galaxies, this is quite high. 
In NGC 660 iBaan et al.l lll992t) find the Hi column density 
Nhi/Ts to be 7.4 x 10^® cm“^/K against their component 
4, the closest to the core, although this is measured with 
the much larger l'.'4 beam of the Very Large Array in A- 
configuration. 


4.2 X-ray spectral fitting 

The low X-ray count rate meant that it was not possible 
to extract a useful spectrum from the core region alone. 
The X-ray spectrum of the entire region of diffuse emission 
around the nucleus was extracted from the data taken after 
outburst (Fig. [HJ. It possesses a continuum which can be 
fairly well fitted by a power law with a spectral index of 
L 6 lloji 2 (reduced — 1.26). 

However, there are several large deviations from this, 
which coincide with well-known emission lines. These can be 
well fitted by Gaussians at 1.80 (^q q 4 keV and 1.35 (!Iq qj keV, 
giving a total reduced of 1-0 together with a power law 
with spectral index 1.59lg j 4 . The widths of these lines, 
5('(g° eV and O.I^q'^j eV respectively, are not well constrained 
due to the low number of counts in the spectra and conse¬ 
quent low spectral resolution. The Gaussians are, however, 
required to account for the excesses in the data and pro¬ 



Figure 6. X-ray spectra from two separate Chandra observa¬ 
tions taken after the radio outburst (red and blue points). The 
line shows a model fitted to the data, consisting of a power law 
and two Gaussians. Both spectra show strong emission lines at 
approximately 1.80 and 1.35 keV, which are most likely due to 
silicon and magnesium transitions. 


vide a good fit. The energies of the Gaussians are consistent 
with emission lines of magnesium and silicon. The features 
are both seen in both post-outburst spectra, and are there¬ 
fore unlikely to be artefacts. These emission lines are com¬ 
monly seen in supernova remnants and are therefore likely 
to be a ssociated with the star for mation observed in this 
region ijGuainazzi fc Bianchil fioo^ ) . Indeed, emission lines 
from these elements have also pr eviously been seen in other 
wavebands llSatvapal et al.ir2004h . 

The Chandra spectra can also be reasonably well-fitted 
with a model for shocked gas which can be used describe 
supernova emission {vpshock, reduced x^ = 1-23), however 
the parameters of the best-fitting models require either a 
high gas temperature (~30 keV) or very high abundances of 
some elements (tens of times higher than solar abundance), 
due to the strength of the lines in the spectra. As would 
be expected, the best fitting parameters require high abun¬ 
dances of silicon and sulphur at low temperatures, as well as 
helium and calcium. As with the images, the pre-outburst 
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Component 

Depth 

Jy 

Width 

kms“^ 

Velocity 

kms~^ 

”^peak 

1 

-0.022 

8.7 

813 

0.32 

2 

-0.016 

2.9 

821 

0.22 

3 

-0.043 

11.7 

838 

0.77 

4 

-0.024 

32.2 

851 

0.36 

5 

-0.024 

3.2 

898 

0.36 

6 

-0.0045 

18.3 

963 

0.06 

1 

-0.017 

10.2 

800 

0.24 

2 

-0.008 

6.3 

816 

0.11 

3 

-0.020 

9.7 

830 

0.29 

4 

-0.008 

4.0 

857 

0.11 

5 

-0.005 

6.8 

894 

0.06 

6 

-0.0007 

18.3 

984 

0.009 


Table 2. Resulting parameters of 6-component fits to the radio 
spectra (Figure |4]l for component 1 (top) and component 2 (bot¬ 
tom). The systematic uncertainty on the velocity measurements 
is ±1.6km/s. 

spectra are of too poor quality to perform meaningful analy¬ 
sis on, but are not inconsistent with possessing strong emis¬ 
sion lines. 


5 THE NATURE OF THE EVENT 

While many radio-loud quasars and galaxies do show some 
degree of time-variable behaviour across the electromagnetic 
spectrum, the ongoing event in NGC 660 is extreme. The 
spectacular and sudden increase in radio emission, together 
with the appearance of new radio structure, has been sug¬ 
gested to originate in either the onset of nuclear activity, 
a tidal disruption event involving the central black hole 
(with precession possibly explaining some of the morphol¬ 
ogy), or a highly luminous supernova event. At Galactic co¬ 
ordinates (141,-47), the likelihood of a foreground Galactic 
flaring event is low, and something this bright, so far from 
the Galactic plane, would be unlikely to have escaped at¬ 
tention. The high-resolution radio structure of this source 
presented here apparently shows a core with two lobes, sim¬ 
ilar to the core-jet morphology often seen in AGN. One likely 
explanation is that the galaxy is undergoing a burst of en¬ 
hanced nuclear activity related to the recent merger or close 
encounter event which created the polar ring and caused the 
starburst activity. In this section we compare the informa¬ 
tion obtained from these new observations and archival data 
with those of possible explanations, some more likely than 
others. 

5.1 A central starburst 

The central compact EVN component has a brightness tem¬ 
perature of Th > 10^ K (see Table [T|. Given that typical 
central starbursts have Tf, ^ 10® K, it is unlikely that an ac¬ 
tive central starburst is a viable explanation for the source. 
Also, since the new source has appeared since the VLBA 
observations in 2001 , and likely since the last reported non¬ 
detection of flaring activity in 2008, the measured brightness 
of the source would imply some 10 to 100 “normal" radio 
supernovae have occurred in a very small region (~l.lpc) 
over a handful of years. This seems unlikely, especially when 


compared to known highly active starbursts. One of the 
most powerful nearby starbursts is that in Arp 220. Here 
[Lonsdale et all (l200a ) find the radio supernova rate to be 
4±2 per year across two nuclei each some 100 to 150 pc 
across, separated by 370 pc. At 13Mpc, 1" is 63 pc, so we 
would easily expect to resolve individual SN in NGG 660 if 
this were the explanation, unless several tens of RSNe oc¬ 
curred within a region of ~lpc. The spectral features de¬ 
tected in the Chandra data do show the presence of chemi¬ 
cals often seen in supernova remnants, although their likely 
presence in the pre-outburst spectra strongly suggest that 
these lines are associated with the pre-outburst star forma¬ 
tion observed in this region, rather than with new activity. 

5.2 A single luminous radio supernova/remnant 

Another possibility is that the outburst is due to a sin¬ 
gle stellar-scale explosive event, although the distance of 
NGC 660 would make this significantly more luminous 
than most known radio supernovae or supernova remnants. 
The 1.4-GHz power from the EVN observation is 1.42 x 
10 ^^ WHz~^; the fact that this is three orders of magnitude 
more powerful than known evolved supernova remnants in 
our o wn Galaxy such as Gas A (10^® WHz“^; iBaars et al.l 
along with the fact that the object waa no more than 
six years old at the time of observation, makes this expla¬ 
nation very unlikely. The highly luminous radio supernovae 
in Arp 220 are a more promising comparison, although even 
here the source in NGC 660 is an order of magnitude more 
powerful ([Lonsdale et al.l [20061 ). 

The extremely energetic GRB-associated supernovae 
SN1998bw and SN2009bb are closer in luminosity to this 
event and resulted in relativistic jets. Both were type Ib/c 
events and relatively short-live d in the radio llKulkarni et al.l 
1 19981 : iBietenholz et al]|2010b[) : the time-scale of the event 
in NGC 660 is more characteristic of a longer-lived type II 
event. One of the brightest RSN ever, SN1986J, has devel¬ 
oped a complex and evol ving morphology, and a spectrum 
which flattens over time (iBietenholz et al]l2010al) . however 
the lower limit on the expansion speed for the source in 
NGC 660 is significantly greater than the expansion seen in 
1986J, or other nearby ordi nary RSN which have been well- 
studied on VLBI scales (e.g. [Bartel et al.ll2002l : [Eenech et al.[ 

[20101 : [Martf-Vidal et alJ[201ll) ' 

A further comparison can be made with the unusual ra¬ 
dio sources in M82. One suggestion for the origin of the fad¬ 
ing source 41.954-575 is a historical GRB seen off-axis from 
our perspective. This source is resolved on VLBA scales and 
would have had a peak luminosity some two orders of mag¬ 
nitude greater than the event in NGC 660, but the VLBI- 
detected component s are expanding at only ~1800km/s 
([MuxIow et al.[[2005[) . 

5.3 A tidal disruption event 

A handful of optically non-active galaxies have displayed 
X-ray flaring b ehaviour interpreted as disruption events 
([Komossa[[2002l ) . The characteristic behaviour of such events 
includes huge peak X-ray luminosity, variability amplitudes 
of up to a factor of ^^200 with a short rise time and declin¬ 
ing over months to years, and an ultra-soft X-ray spectrum. 
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In the case of NGC 660, the expected X-ray peak luminos¬ 
ity should have exceeded the flux of the diffuse continuum 
by several orders of magnitude, but the archival X-ray data 
does not show this expected huge peak in luminosity. 

An example of an X-ray flare thought to be due to 
a stellar disruption event is the case of NGC5905. An 
X-ray flare was detected in 1991 and VLA observations 
were carried out to search for post-outburst radio emis¬ 
sion in 1996, but no compact central object was detected 
putting a limit of < 10^° WHz~^ on the radio luminosity 
jKomossa fc Dahlemll200lh . On e-MERLIN scales we see a 
source appear in NGC 660 compared to archival data at the 
same resolution, implying a variability factor of ~1200 in 
the radio. However the rise time is unknown, and the decay 
is, as yet, not well-sampled. There is no evidence for prior 
low-luminosity AGN activity, with no optically-hidden radio 
AGN present in archival radio data. 

The case of Swift J1644+57 also provides an interest¬ 
ing comparison. Displaying a mildly relativistic outflow and 
spectral f eatures due t o sync hrotron and inverse Compton 
processes, iBloom et al.l (1201 il l draw an analogy with a tem¬ 
porary small-scale blazar-like event. iGiannios fc Metzgeil 
(1201111 discuss transient radio emission from such events 
and model the resulting characteristics, focussing on off-axis 
events. In the case of NGC 660, we have no information on 
the time-scale from outburst to peak luminosity, and infor¬ 
mation on the rate of decay will only come with time. 


5.4 A new period of AGN activity 


Whilst low-power jets from tidal disruption events are short¬ 
lived and dissipate their energy in to the ISM rapidly, jets 
from AGN last orders of magnitude longer and propagate 
to much larger distances. The most likely scenario is that 
the central super-massive black hole is beginning a period 
of activity, turning NGC 660 from a LINER-type galaxy into 
a low-luminosity AGN. 

There is now much evidence that a significant frac¬ 
tion of nearby LLAGN galaxies host flat-spectrum compact 
cores with associated parsec-scale jets ( e.g. in M51 Ram- 
pada r ath et al in prep; more generally in Nagax_eLal 20021 . 
I2OO5I : IPanessa fc Girolettil 12013 ; iMezcua fc Prietdl2014l L~in 
NGC 660, the central source detected by e-MERLIN is both 
compact and flat spectrum, and the EVN map shows ev¬ 
idence of a core-jet structure. The spectral index of the 
EVN sources cannot be calculated in this case, since there 
is only 1.4-GHz information available at present. The EVN 
morphology resembles that of a compact symmetric object. 
These are the youngest class of AGN, with steep spec¬ 
tra, measurable expansion and the lobes entirely contained 
within the (clumpy, asymmetric) medium of the host galaxy. 
Such activity is generally triggered by large-scale collision or 
merger events driving material into the central region of the 
galaxy, and NGG660 has undergone a relatively recent en¬ 
counter which greatly disturbed its morphology. 

Whilst limits placed on the apparent expansion velocity 
by the archival data suggest that the motion could be mildly 
relativistic, if components 2 and 3 are a jet and counter-jet 
respectively, then the ratio of their flux densities (-^2.5) sug¬ 
gests that they are not moving relativistically. Eurther EVN 
observations have been requested in order to look for mor¬ 
phological changes or relative motion of the components. If 


1.5 and 5 GHz e-MERLIN May/June 2013 



Figure 7. Peak flux measurements per sub-band from e-MERLIN 
observations made at 1.5 and 5 GHz in May/June 2013 (-I-). The 
lines shows least squares fits to the e-MERLIN data, with solid 
black, dashed black and dotted grey lines corresponding to models 
(i), (ii) and (iii) in the text (see Section l5.4ll . 


this source is a newly-formed jet, these observations will en¬ 
able us to better constrain the expansion velocity, resulting 
in a better estimate for the onset of activity. 

As discussed in section [S] the average spectral index of 
the new radio component is a = 0.5. With the wide band¬ 
width of the new e-MERLIN system it is possible to measure 
the spectral index within each band; EigureQshows the flux 
measurements for the new (unresolved on e-MERLIN scales) 
component in each of the eight sub-bands at 1.4 GHz and 
four sub-bands at 5 GHz. The results clearly show that this 
new radio object is a GHz-peaked spectrum (GPS) source, 
with a synchrotron spectrum above a turnover at around 
3—3.5 GHz. The spectral index below the peak is a = -1-1.6, 
while above the pe ak it is a = —0.72. Erom a large sam¬ 
ple of GPS sources, Ide Vries et al.l lll997ll found that (using 
the same convention S' oc j/“) the average spectral index on 
the optically thick side of the turnover frequency is 4-0.51, 
whilst that on the optically thin side is —0.73. NGC 660 is 
steeper on the optically thick side than this average, suggest¬ 
ing stronger self- absorption, although there is a large scatter 
in the sample of ide Vries et al.l (ll997l L 

GPS sources generally represent a very early phase in 
the evolution of radio galaxies, with their extent entirely 
contained within the central kpc of their host galaxy, mor¬ 
phologies similar to those of larger more evolved classical 
double-lobed radio galaxies, and the pea k in their spe ctra 
likely due to synchrotron self-absorption (IO’Deal[l998l and 
refs therein). This description closely matches the observed 
properties of the ~6-year old source in NGC 660. A young 
source would be expected to have a steep, constant spec¬ 
tral index, while older AGN-related structures (jets, lobes) 
flatten at low frequency due to radiative losses. In the case 
of NGC 660, the EVN observations clearly show multiple 
source components within the e-MERLIN beam. 

Fitting the e-MERLIN spectral points with (i) a sim¬ 
ple power law, (ii) a power law with free-free absorption by 
a foreground screen of ionised gas, and (iii) a self-absorbed 
Bremsstrahlung spectrum, we find that none of these models 
are a good fit to the data, all being less steep at frequen¬ 
cies above the turnover (Fig. [7|). Further e-MERLIN mea- 





































8 M.K. Argo et al 


surements have been requested, covering a wider frequency 
range at C-band (4—8GHz), to better constrain the fits, but 
the explanation could simply be that the different compact 
components seen with the EVN have differ ent spectra, an 
effect also discussed bv ide Vries et all lll997ll in the context 
of compact double sources. 

The X-ray spectra of the emission region after the out¬ 
burst show a total flux of 1.85^g x 10~^^ er gscm~ ^ s~^. 
Calculations using the method from iMarsch^ Jl983ll show 
that, if the radio emission seen at the nucleus is due to syn¬ 
chrotron self-absorption, the X-ray flux expected from the 
region is negligibly small for an assumed radio nucleus size 
of ~0.5mas, from the VLBA observations. All of the pa¬ 
rameters in this calculation are either well constrained or 
do not have a significant impact on the final value, except 
for the source size. Since the core itself is unresolved, the 
physical size could be smaller. A typical core size of 0.2 mas 
would still produce only very low X-ray fluxes. A size of 
~0.1 mas or s maller, which is the case in only a small frac¬ 
tion of AGN llHoriuchi et al.ir2004h . is needed to produce a 
flux comparable to the 10“^® erg cm“^ s”*^ diffuse flux. If the 
source is a newly-awoken AGN, X-ray emission from addi¬ 
tional processes to synchrotron self-Gompton are likely to 
be present, however, many examples of low luminosity AGN 
with X-ray luminosities of < ergs“^ exist, equating to 
a flux of < 10~^^ ergc m~^ s~^ at the distance of NGC 660 
llTerashima &: WilsorJ 120031 . Furthermore, the radio spec¬ 
trum indicates a very high absorbing column in the nucleus, 
which would cause a large degree of X-ray absorption if cov¬ 
ering the X-ray source. It is therefore not unexpected that 
no strong X-ray nuclear emission is observed, even if a new 
AGN has appeared. If the increase in flux is due to a new 
period of SMBH activity in the centre of NGC 660, then the 
radio luminosity measu red here places it re latively low on 
the fundamental plane dMerloni et al.ll2OT3l l. hence the X- 
ray luminosity would also be expected to be low and the 
non-detection is consistent. 


6 CONCLUSIONS 

The outburst in NGC 660 has been explored using archival 
radio and X-ray data, together with new observations at 1.4 
and 5 GHz with e-MERLIN, and at high resolution with the 
EVN at 1.4 GHz. Our observations show the appearance of 
a bright new source in the core of the galaxy, an object not 
present in any archival radio images of the galaxy prior to 
2010. This new object is compact on a scale of ~50 milliarc- 
seconds with e-MERLIN, but displays an obvious jet-like 
component to the north-east and a weaker possible counter¬ 
jet to the west at higher angular resolution with the EVN. 
The evidence from these new radio observations, combined 
with archival data from Chandra, MERLIN, the VLA and 
the VLBA, suggests that we are seeing the activation of a 
low-luminosity AGN in NGC 660. 

Hi spectra obtained during the EVN observation show 
strong absorption along this new line of sight through the 
disturbed disk of the galaxy, potentially absorbing any asso¬ 
ciated X-ray emission. The Hi spectrum of the new central 
radio source strongly resembles that seen at lower resolution 
with Westerbork-only data, post-outburst, with several nar¬ 
row velocity components due to absorption in the disk, and 


a slightly broader red-shifted wing suggesting an inflow of 
material along our line of sight. 

The EVN flux point at 1.4 GHz is almost 50% lower 
than that measured by e-MERLIN at the same frequency 
130 days earlier. This is likely due to a combination of over¬ 
all fading of the source, spectral evolution, and resolution 
effects between the two arrays; the relative effects of these 
two explanations will be tested by planned follow-up ob¬ 
servations. The astrometric capability of the EVN will also 
allow us to better constrain the direction of expansion and 
the velocity of the components, resulting in a better esti¬ 
mate for to; an apparent velocity of 1.2c would result in a 
movement of 0.4 pc over one Earth-year, a shift of 5.8 mas 
at a distance of 13 Mpc. There is also much information to 
be extracted from further e-MERLIN observations at lower 
angular resolution. Observations at both 1.4 and 5 GHz have 
been requested in order to measure any changes to the ra¬ 
dio spectral energy distribution of the new source, and ad¬ 
ditionally probe several additional spectral lines, providing 
valuable Information on physical conditions of the gas in the 
centre of NGG 660. 

Whatever the cause of the increase in radio luminos¬ 
ity, the fact that such variability does occur in LINER-type 
galaxies is potentially a serious issue for statistical studies 
of the link between star formation, fuelling and accretion in 
low luminosity AGN. If this is, as it appears to be, a new 
period of AGN activity, then in the short term we can ex¬ 
pect to track the expansion as the jets move out through the 
galactic medium, and watch for evolution of the GPS-nature 
of the spectrum. In the much longer term, if strong enough 
the outflows could quench future star formation by clearing 
out the gas reservoirs needed. GPS sources are seen across 
a wide variety of peak luminosities, galaxy types and red- 
shifts, and are useful probes of the early universe. NGC 660 
is now one of the closest examples, and presents an ideal 
opportunity to watch as it evolves. 
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